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EXACT SOLUTIONS AND PARTICLE CREATION FOR 

NONCONFORMAL SCALAR FIELDS IN HOMOGENEOUS 

ISOTROPIC COSMOLOGICAL MODELS 

Yu. V. Pavlov 

"en ' 

The problem is solved of describing scale factors of a homogeneous isotropic 
£NJ ' spaces-time such that the exact solution for the scalar field with a nonconformal 

coupling to curvature can be obtained from solutions for the conformally coupled 
field by redefining the mass and momentum. Explicit expressions for dependence 
of time from the large-scale factor are presented in the form of Abelian integrals in 
these cases. The exact solution for a scalar field with Gauss-Bonnet type coupling 
with curvature is received and it is shown that the corresponding nonconformal 
additions can dominate at the particles creation by gravitational field. 
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^ ! 1. Introduction 

oo 



Quantum effects in a curved space-time, in particular, particle creation by the grav- 
itational field, can have important applications in cosmology and astrophysics PQ, [2]. 
Quantum theory in a gravitational field has been studied intensely since the 1970s. The 
case studied in particular detail is that of conformally coupled fields in a homogeneous 
^ isotropic space-time. Less studied is the case of a nonconformal coupling of a scalar field to 

curvature. At the same time, nonconformal contributions can dominate both in the effect 
of particle creation and in the magnitude of vacuum averages of the energy-momentum 
tensor, which determinate the back reaction of a quantized field on the space-time metric 
(see, e.g., [3], [1]). Until recently, no calculation of quantum effects of a scalar field coupled 
to a curvature of the Gauss-Bonnet (GB) type have been available. 

Here, we study the possibility of obtaining exact solutions for scalar fields with both 
the standard nonconformal and the GB couplings to curvature from solutions in the case 
of a conformal coupling in a homogeneous isotropic space-time. We review scaling factors 
that allow deriving exact solutions for nonconformal fields. We give an exact solution for 
the scalar field coupled to the GB-type curvature and analyze particle creation in such a 
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model. We calculate particle creation in the quantum field theory in a curved space-time, 
where quantized matter fields are considered in an external classical gravitational field 
described by the space-time metric. The possibility of using such a semiclassical approach 
relates to a broad domain between Planckian and Compton characteristic values of length, 
curvature, and density p], [2]. 

We use the system of units where h — c—1. 

2. Scalar field in a homogeneous isotropic space 

We consider a complex scalar field if (x) of mass m with the Lagrangian 

L(x) = J\g~\ [g ik d i( p*d k( p - (m 2 + V g ) <p*cp] (1) 

and the corresponding equation of motion 

(V 4 V, + V g + m 2 ) ip(x) = , (2) 

where V» are covariant derivatives in the N- dimensional space-time with the metric g^, 
g = det(gik), and V g is a function of invariant combinations of the metric tensor g^ 
and its partial derivatives. The case V g = corresponds to the minimal coupling of 
the scalar field to curvature. The case V g = £ C R, where R is the scalar curvature and 
£ c = (N — 2)/ [4(iV — 1)], ia called the conformal coupling to curvature (£ c = 1/6 at 
N — 4). Equation (J5J) is conformally invariant if m = and V g = £ C R. 

Taking an arbitrary V g leads to the appearance of third- and higher-order derivatives 
of the metric in the metric energy-momentum tensor of the scalar field and hence in the 
Einstein equations. The appearance of additional higher-derivative terms in the equations, 
even if their coefficients are small, leads to a dramatic restructuring of the theory. If we 
require that the metric energy-momentum tensor not contain derivatives of the metric of 
order higher than two, then V g can be chosen as the function 

V g = ^R + CRcB, (3) 

where 

p 2 del T-) T-) Impq j\t) T)lm , p2 

K GB — K lmpqti y — ^rilraJX + K . 

In the four- dimensional space-time, Rq B coincides with the Euler characteristic density, 
which by the GB theorem is a topological invariant (in N = 2, the corresponding density 



is proportional to R). Therefore, the coupling to a curvature of form ([3]), introduced 
in [5], can be naturally called the GB-type coupling. In lower dimensions N = 2,3, 
new effects due to ( ^ are absent, because then Rq B = 0. For N = 4 with constant 
(f(x), the contribution of the Rqb term to the metric energy- momentum tensor is absent 
because the corresponding variational derivative vanishes [6j. But for a variable <p(x), the 
contribution od such terms must be taken into account if the constant ( of dimension 
(mass) -2 is nonzero. 

We write the metric of an iV-dimensional homogeneous isotropic space-time in the 
form 

ds 2 = dt 2 - a 2 (t) dl 2 = a 2 (r]) (drf - dl 2 ) , (4) 

where dl 2 = / j a/ 3dx a dx^ is the metric of the (iV— l)-dimensional space of constant curvature 
K = 0, ±1. The full system of solutions of Eq. (T5]) in metric @ can be found in the form 

^ x )=a^ N - 2 ^ 2 {r 1 )gM^J^), (5) 

where 

g'i( V ) + n 2 (r ] )g x (r ] ) = 0, (6) 

n 2 (n) =m 2 a 2 + \ 2 -A£a 2 R + (a 2 R^ Bj (7) 

A^^x) = -(a 2 - (^) Yj-Mx), (8) 

A£ = £ c — £, the prime denotes the derivative with respect to the conformal time r], and J 
is the set of indices (quantum numbers) that label eigenfunctions of the Laplace-Beltrami 
operator Aat_i in the (N — l)-dimensional space. In metric (Hj), the expressions for the 
scalar curvature and the GB invariant are [S] 

R = a - 2 (N-l)[2c' + (N -2)(c 2 + K)], (9) 

R* B = a-\N-l)(N-2)(N-3)(c 2 + K) [Ac' + (iV-4) (c 2 + K) ] , (10) 

where c = a' /a = d(t). 

In accordance with the method of diagonalization of a Hamiltonian pQ (see [7] for the 
case of an arbitrary V g ), the functions gx(v) must satisfy the initial conditions 

g' x (vo) = itt(vo)g\(vo), \g\(vo)\ = ^ 1/2 (vo)- ( n ) 



If the quantized scalar field is in the vacuum state at the instant 770, then the number 
density of the pairs of particles created up to the instant r\ can be evaluated (for K — 0) 
as [1] 



X 



»(»?) = 2^T ISM* N -*d\, (12) 



where B 



N 




-1 



2 iv-3 7r (Ar-i)/2 F (( A r_ x y 2 ) , T(z) is the gamma function, and 

SM= Mn)-in 9 M\ 2 (13) 

As shown in [8], S*a ~ A -6 , and the integral in ( TT2l) converges for N < 7. 



3. On exact solutions for different types of coupling to 

curvature 

With the substitution g{rj) = expz(r/), Eq. ([6]) reduces to the Riccati equation of the 
general type for v{rj) = z'{rf): 

v'(7 1 )+v 2 (r ] ) + n 2 (7 1 ) = J (14) 

Therefore, the number of scaling factors a(rj) admitting exact solutions is relatively small. 
In the cases where an exact solution can nevertheless be found, it is typically expressed 
in terms of special functions: hypergeometric, Bessel functions, etc. Several exact solu- 
tions are given in [TJ, [2]. Scaling factors admitting exact solutions are briefly reviewed 
in 0, [ID]. 

In a homogeneous isotropic space-time, Eq. ([6]) has a simpler form for the field with 
the conformal coupling Q 2 = m 2 a 2 + A 2 than for a nonconformal coupling. This does not 
forbid the existence of exact analytic solutions of a more complicated equation for a non- 
conformal scalar field. But the majority of the known exact solutions for a nonconformal 
scalar field are obtained from solutions for the conformal field by replacing the mass and 
momentum parameters. 

We consider the following question. At what scaling factors can the exact solution for 
the scalar field with a nonconformal coupling to curvature be obtained from a solution 
for the field with conformal coupling by redefining the values of the field mass m and the 
dimensionless momentum A? 



For the scalar field with a nonconformal coupling of the form £R, in accordance 
with (J7J), this is to be the case if there exist constants /3 X and (3 2 suc h that 

A£Ra 2 =p 1 +p 2 a 2 . (15) 

This condition, with fl9]) taken into account, is an ordinary second-order differential equa- 
tion for a(rj), which reduces to a linear first-order equation for f(a) = a' 2 (r]). Its solution 
can be represented in the form 

da 



/ r , ,_ v r , 9^^ 4 =± (v~Vo), (16) 

VCoo 4 JV + Cia z + 6 2 a 4 

where Co, C\, C 2 , and tjq are arbitrary real constants such that C$a A ~ N + C\a 2 + C*2a 4 > 0. 

In the general case for N > 4, integral ( 1T6J) is hyperelliptic, and calculating 0(77) 
amounts to inverting the hyperelliptic integral [II] . For iV = 4, obtaining 0(77) from ( fT6l) 
is the problem of inversion an elliptic integral p2] . In particular cases where integral ( TT6i) 
is pseudoelliptic or N = 4 and C2 = 0, it can be expressed in elementary functions. We 
give examples of such scaling factors. 

The first example of a scaling factor is given by the relations 

C 1 = C 2 = 0, a = ai r] 2/{N - 2) = a t 2/N , a 2 R = (N - 1)(N - 2)K, (17) 

where a\, 00 = const. For N = A with K = 0, such a scaling factor corresponds to a case 
extremely important from the application standpoint, the radiation-dominated universe. 
Solutions are known for N = 4, where they can be expressed in terms of Kummer's 
hypergeometric function, and for N = 6, where they can be expressed in terms of Bessel 
functions [9], [T3] . 

The second example of a scaling factor is 

C = 6*2 = 0, a = ai e ari = at, a 2 R = (N - 1)(N - 2) (a 2 + K), (18) 

where a = const. For a = 1,K = —1, this is Milne's universe [14] (i.e., a part of 
Minkowski space in the appropriate coordinates). A solution of Eq. (EJ) expressed in 
terms of a Bessel function was presented in [5]. 
The third example of a scaling factor is 

-1 e Ht r 1 

C = d = 0, C 2 = H 2 , a = _ = _ a 2 R=(N-l) NH 2 a 2 + (N-2)K , (19) 

Hr] H L J 



where H = const. For K = 0, corresponding coordinates describe a part of the De Sitter 
universe. The solution of Eq. fl6]) can be expressed in terms of a Hankel function [TJ. 
For C = 0, solutions of Eq. (JTBjl are also given by 

01 . it -oi t* ai 7t 

ai sm — , -— = a\ sinn — , = oi cosh — , (20) 



cosh 777 a\ sinh7r/ a\ cos 777 a\ 

where 7 = const. In the four-dimensional space-time, solutions of Eq. (1161) with C\ 
4CqC2 are 



2jt / -27t 
ai tan 777 = aiWexp 1, ai tanh7?7 = an/1 — exp . 

V a i V a i 



27* < 

— a\ coth7r7 = aiWexp h 1. 

V ai 



(21) 



and solutions with C<i = are 



ai sin7?7 = yaf — 7 2 t 2 , ai sinh7^ = y7 2 t 2 — a 2 , ai cosh 777 = y7 2 t 2 + a 2 , (22) 

Solutions of Eq. ([6]) with scaling factors (120]) and ff2T]) can be expressed in terms of the 
hypergeometric functions. Equation ([6]) with scaling factor ( )22|) reduces to a (modified) 
Mathieu equation. 

We note that if a solution of Eq. ([6]) with a scaling factor a{rf) is found, then it can 



be used to obtain a solution for the scaling factor a = J a?{rj) + b 2 , where b = const, by 
redefining A. 

With the coordinate time t, solutions of Eq. ffT5|) can be represented in the form 



/ 



d (a N ' 2 ) N 

-±-(t-t ), (23) 



Co + C iC l n - 2 + C 2 a N 2 



where to — const. For N = 4, integral ( 1231) is expressed in elementary functions. There- 
fore, in the four- dimensional space-time, the dependence of the coordinate timet on scaling 
factor under which solutions for the scalar field with a £i? coupling to curvature can be 
obtained from solutions for the field with conformal coupling by redefining the mass and 
momentum values must be expressible in a finite form in terms of elementary functions. 
We consider the case of a scalar field with the GB coupling to the curvature in (^ and 
assume that N > 4 and ( 7^ 0. An exact solution of Eq. ([6]) can be obtained from the 



solution for the conformally coupled field by redefining m and A if there exist constants /3i 
and 02 suc h that 

- A£R a 2 + (Rl B a 2 = fa + /3 2 a 2 . (24) 

With formulas fl9]) and ffTUj) . this condition reduces to a first-order linear inhomogeneous 
differential equation for the function 

na)=^(t) + K- 2(iN A lf {N _ 3) ) 2 . (25) 

Solving differential equation (j2"4"|) for a(t), we obtain an expression for the dependence of 
the scaling factor on the coordinate time: 

' ia = ±(t - to), (26) 



where D , Di, and D 2 are arbitrary real constants such that the integral is real. 

In the general case, finding the explicit form of a(t) amounts to inverting Abelian 
integral ( 1261) . We give examples in particular cases: 

K = A£ = D 1 = D 2 = 0, a = 7 t 4/7V , R% B = 0, (27) 

where exact solutions for the conformal coupling are known for N — 4, 8, 12 (see formu- 
las ([I7D and pgjl ); 

AT = D = D 1 = 0, a = a expHt } R 2 B = N(N-1)(N-2)(N-3)H 4 (28) 

(see formulas (fl9l) ); and for N — 4, functions (!20l) are solutions of Eq. (|26l) with D 2 = 
AD Q D 2 . 

In what follows, we consider the exact solution and particle creation for the first scaling 
factor in (J20j). 

4. Particle creation in the a\ smi^t/di) = a\j cosh 77/ model 

Space-time with a scaling factor a\ sin(7t/ai) = ai/ cosh 777 evolves between two sin- 
gularities at t — (77 = —00) and t = T = nai/j (77 = +00), if 7 2 + if 7^ 0, iV > 3. In 
the case 7 = 1 and K = —1, the corresponding coordinates cover a part of the De Sitter 
space of the second kind [15] with a constant negative curvature R — —(N — l)N/a 2 . 
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Particle creation in this model was previously considered for a conformal coupling 
to the curvature in pQ, and the case of the £,R coupling was investigated in [3]. We here 
consider the creation of scalar particles with a GB-type coupling to the curvature, Eq. (|3J), 
in a four-dimensional space-time. From fl9]) and ffTOj) with N = 4, we obtain 



Ra 2 = 6 ( 7 2 + K 



12 7 2 



a K GB 



24 7 2 



^n 1 -^ 



and therefore Eq. ([6]) becomes 



</"fa) + 



■i > aho 2 , ,24 7 4 ^ 



m^+Aei2y + C- 



fl i 



' + A 2 - 6( 7 2 + if) (A£ + C 47 '^ 



(29) 

0(77) = 0. 
(30) 
The solution of Eq. (I3"0~|) with initial conditions (ITTj) with 77 — >• — 00 can be obtained from 
the exact solution for the conformally coupled case pQ by redefining the values of the mass 
(m — ¥ M) and momentum (A — >• A): 



cosh 777 



a\ 



M 1 






«I 



A 2 = A 2 -6( 7 2 + if)(A£ + C^ 



The solution has the form 

g(v) 



o i(Ar)+a ) 



F[A,B;C; 



1 + tanh 777 



(31) 



(32) 



'A V ' ' ' 2 

where -F(A, B; C; z) is the hypergeometric function, «q is an arbitrary real constant, and 



.4 



1 



1 M 2 a\ 



B 



\4 



1 M 2 a\ 



C = 1 + i 



A 



T 



7 



(33) 



2 \| 4 7 

The condition f2 2 > 0, which is necessary in the Hamiltonian diagonalization method 
is satisfied for Eq. ( 130]) for any A and rj if 



(7 2 + if) A£ + C 



4T 2 
a 2 



< 0, m 2 a\ + 6(7 2 - K)A£ - 24(k\ > 0. 



In particular, Q 2 > 0, at if = and 



4 7 2 
of 



2^2 



ma 



6 7 2 



(34) 



(35) 



From asymptotic expression of ({TBI for exact solution (1321) . similarly to the conformal 
case pp, we obtain the limit spectrum (as 77 — > +00) of the created particles: 



71 



COS 



s> 



W' a 

sinh — 

7 



4M 2 a 2 



V 



(36) 



In accordance with formulas ( TT2|) . (l3Tj) . and (1361) . the number of particle pairs created in 
a Lagrangian volume a 3 (t) at K = per the evolution cycle t — > T is given by 



AT 



7 2 

COS 71% 



2tt 5 



2 °r X\/x 2 — Xn 



\ 4 7^ af i smh x 

N x 



<ic, 



where 



Because [TO 



x = irJ -6 (&£ + (% 



X 2 IT 2 

sinh x 6 



it follows from formulas ( 131"]) . ( |35l) . and ( )38|) that 



(37) 



(3? 



(39) 



>47" ^ m"a 



2„2 
1 



a l 



67 s 



AT 



7 2 71 " 

12^ COS 2^ 



1 + 24A£ + 



4m 2 a 2 

ry2 



(40) 



Formulas (IBT)]) and ( 13TJ) reproduce the corresponding results for the scalar field with 
the £R coupling at ( = [3] and for conformally coupled particles at ( = = A£ pQ. At 
A£ = 0, formula ( HD1 gives the number of pairs of conformally coupled particles created 
per the evolution cycle: 



Nn 



T 



cos" — < 



1 + 



4m 2 a 2 



(41) 



12tt 3 "'" 2^ Y 

If ( — y — oo for a fixed A£, then M 2 — > — oo, and the number of created particles 
increases indefinitely, as is shown, for example, for moi/7 = 1 in Fig. [TJ 
It follows from fl36l) that particle creation per evolution cycle is absent if 



M l a{ = n(n + 1)Y 



n = 0,1,2,... 



(42) 



With expression f l3~T|) and inequalities fl35|) with AT = taken into account, there is no 
particle creation per evolution cycle if the parameter ( takes the values 



c. 



24 7 2 



n{n + 1) - 12A£ 



m 2 a 2 i 



n = 0, !,...< 



m 2 a 2 



N ^ + 6Af + i-i. (43) 



Hence, depending on the values of (, particle creation in this model can either increase 
without bound or be absent altogether. Therefore, the effect of the parameter ( of the 
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1000 




tf 



Figure 1: Dependence of the number of created particles on the parameter (: 
curve 1, A£ = —1/6; curve 2, A£ = 1; and curve 3, A£ = 3. 

scalar field coupling to the GB invariant can dominate, and it must be taken into account 
in calculating of the effects of scalar fields in a curved space-time. 
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